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Abstract 

The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic de- 
scription of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field 
crystal model describes a variety of spatially localized structures, in addition to different spatially extended 
periodic structures. The location of these structures in the temperature versus mean order parameter plane 
is determined using a combination of numerical continuation in one dimension and direct numerical simu- 
lation in two and three dimensions. Localized states are found in the region of thermodynamic coexistence 
between the homogeneous and structured phases, and may lie outside of the binodal for these states. The 
results are related to the phenomenon of slanted snaking but take the form of standard homoclinic snaking 
when the mean order parameter is plotted as a function of the chemical potential, and are expected to carry 
over to related models with a conserved order parameter. 
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I. INTRODUCTION 



Spatially localized structures (hereafter LS) are observed in a great variety of pattern-forming 
systems. Such states include spot-like structures found in reaction-diffusion systems [1] and in 
a hquid light-valve experiment [2], isolated spikes observed in a ferrofluid experiment [3] and 
localized buckled states resulting from the buckling of slender structures under compression [4]. 
Related states are observed in fluid mechanics, including convection [5-7] and shear flows [8]. In 
other systems, such as Faraday waves, the LS oscillate in time, either periodically or with a more 
complex time-dependence, forming structures referred to as oscillons [9, 10]. This is also the case 
for oscillons in granular media [11] and in optics [12]. Other examples of LS include localized 
traveling waves [13, 14] and states called "worms" observed in electroconvection [15, 16]. In many 
of these systems the use of envelope equations removes the (fast) time-dependence and maps such 
time-dependent structures onto equilibria. 

Many of the structures mentioned above are examples of "dissipative solitons" [17] in which en- 
ergy loss through dissipation is balanced by energy input through spatially homogeneous forcing. 
Others (e.g., [3, 4]) correspond to local minima of an underlying energy or Lyapunov functional. 
This is the case for systems whose dynamics are of gradient type, provided only that the energy 
is bounded from below. On finite domains with null boundary conditions such systems evolve to 
steady states which may or may not correspond to LS. However, in generic systems supporting 
dissipative solitons no Lyapunov function will be present and the evolution of the system may be 
nonmonotonic. 

In either system type, time-independent LS are found in regions of parameter space in which 
a time-independent spatially homogeneous (i.e., uniform) state coexists with a time-independent 
spatially periodic state. Within this region there is a subregion referred to as the snaking or pinning 
region [18] in which a great variety of stationary LS are present. In the simplest, one-dimensional, 
case the LS consist of a segment of the periodic state embedded in a homogeneous background; 
segments of the homogeneous state embedded in periodic background may be thought of as local- 
ized hole states LH. Steady states of this type lie on a pair of intertwined branches that snake back 
and forth across the snaking region, one of which consists of reflection-symmetric LS with a peak 
in the center, while the other consists of similar states but with a dip in the center, and likewise 
for the holes. Near the left edge of the snaking region each LS adds a pair of new cells, one at 
either end, and these grow to full strength as one follows the LS across the snaking region to its 

2 

Prenrint- contact: u.thiele@lboro.ac.uk - www.uwethiele.de - January 21 . 201 3 



right edge, where both branches turn around and the process repeats. Thus, as one follows the LS 
branches towards larger norm, both types of LS gradually grow in length, and all such struc- 
tures coexist within the snaking region. On a finite interval the long LS take the form of holes in an 
otherwise spatially periodic state but on the real line, the LS and LH remain distinct although both 
occupy the same snaking region. The LS branches are, in adition, interconnected by cross-links 
resembling the rungs of a ladder, consisting of asymmetric LS [19]. In generic systems posed on 
the real line states of this type drift, either to the left or the right, depending on the asymmetry, but 
in systems with gradient dynamics the asymmetric states are also time-independent. These, along 
with bound states of two, three, etc LS/LH are also present within the snaking region [19]. 

The above behavior is typical of nonconserved order parameter fields. However, an important 
subclass of gradient systems possesses a conserved quantity, and in such systems the order pa- 
rameter field has a fixed mean value. Systems of this type arise frequently in fluid convection and 
other applications [20-22] and are distinguished from the standard scenario summarised above 
by the following properties [23, 24]: (i) the snaking becomes slanted (sometimes referred to as 
"sidewinding"), (ii) LS may be present outside of the region of coexistence of the homogeneous 
and periodic states, (iii) LS are present even when the periodic states bifurcate supercritically, i.e., 
when the coexistence region is absent entirely. The slanting of the snakes-and-ladders structure 
is a finite size effect: in a finite domain expulsion of the conserved quantity from the LS implies 
its shortage outside, a fact that progressively delays (to stronger forcing) the formation events 
whereby the LS grow in length. The net effect is that LS are found in a much broader region of 
parameter space than in nonconserved systems. 

The above properties are shared by many of the models arising in dynamical density functional 
theory (DDFT) and related phase field models of crystalline solids. The simplest such phase field 
crystal (PFC) model [25] (see below) leads to the so-called conserved Swift-Hohenberg (cSH) 
equation. This equation was first derived, to the authors' knowledge, as the equation governing 
the evolution of binary fluid convection between thermally insulating boundary conditions [26]; 
for recent derivations in the PFC context see Refs. [25, 27, 28]. In this connection the PFC model 
may be viewed as probably the simplest microscopic model for the freezing transition that can be 
constructed. In this model the transition from a homogeneous state to a periodic state corresponds 
to the transition from a uniform density liquid to a periodic crystalline solid. The LS of interest 
in this model then correspond to states in which a finite size portion of the periodic crystalline 
phase coexists with the uniform density liquid phase, and these are expected to be present in the 
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coexistence region between the two phases. Some rather striking examples of LS in large two- 
dimensional systems include snow-flake-like and dendritic structures, e.g., Refs. [25, 29-32]. In 
fact, as shown below, the LS are also present at state points outside of the coexistence region. 
However, despite the application of the cSH (or PFC) equation in this and other areas, the detailed 
properties of the LS described by this equation have not been investigated. In this paper we make 
a detailed study of the properties of this equation in one spatial dimension with the aim of setting 
this equation firmly in the body of literature dealing with spatially localized structures. Our results 
are therefore interpreted within both languages, in an attempt to make existing understanding of 
LS accessible to those working on nonequilibrium models of solids, and to use the simplest PFC 
model to exemplify the theory. In addition, motivated by Refs. [29-32], we also describe related 
results in two (2d) and three (3d) dimensions, where (many) more types of regular patterns are 
present and hence many more types of LS. Our work focuses on 'bump' states (also referred to as 
'spots') which are readily found in direct numerical simulations of the conserved Swift-Hohenberg 
equation as well as in other systems [2]. 

Although the theory for these cases in 2d and 3d is less well developed [33, 34] continuation re- 
sults indicate that some of the various different types of LS can have quite different properties. For 
example, the bump states differ from the target-like LS formed from the stripe state that can also be 
seen in the model. In particular, spots in the nonconserved Swift-Hohenberg equation in the plane 
bifurcate from the homogeneous state regardless of whether stripes are subcritical or supercriti- 
cal [35], see also Ref. [36]. The key question, hitherto unanswered, is whether two-dimensional 
structures in the plane snake indefinitely and likewise for three-dimensional structures. 

The paper is organized as follows. In Sec. II we describe the conserved Swift-Hohenberg 
equation and its basic properties. In Sec. Ill we describe the properties of LS in one spatial 
dimension as determined by numerical continuation. In Sec. IV we describe related results in 
two and three spatial dimensions, but obtained by direct numerical simulation of the PFC model. 
Since this model has a gradient structure, on a finite domain all solutions necessarily approach a 
time-independent equilibrium. However, as we shall see, the number of competing equilibria may 
be very large and different equilibria are reached depending on the initial conditions employed. In 
Sec. V we put our results into context and present brief conclusions. 
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II. THE CONSERVED SWIFT-HOHENBERG EQUATION 



A. Equation and its variants 

We write the cSH (or PFC) equation in the form 

9,0(x,t) = «V^-^, (1) 

where 0(x, t) is an order parameter field that corresponds in the PFC context to a scaled density 
profile, a is a (constant) mobility coefficient and F[(f)] denotes the free energy functional 



= J 



(2) 



Here x = {x,y,z), V = {dx,dy,dz)^ is the gradient operator, and subscripts denote partial 
derivatives. It follows that the system evolves according to the cSH equation 

dt(p = aV' [r0 + (V' + + 0^] . (3) 

Equation (3) is sometimes called the derivative Swift-Hohenberg equation [37, 38]; many pa- 
pers use a different sign convention for the parameter r (e.g., [27, 29, 39-45]). In this equation the 
quartic term in F[(f)] may be replaced by other types of nonlinearity, such as = — 620'^/3 + (p^/ 4 
[27, 29, 37, 38] without substantial change in behavior. Related but nonconserved equations 
dt(f) = —aSF[(j)]/6(j) with nonlinear terms of the form /23 [46] or /as = —630^/4 + 0^/6 [19] 
have also been extensively studied, subject to the conditions 62 > 27/38 (resp., 63 > 0) required 
to guarantee the presence of an interval of coexistence between the homogeneous state = and 
a spatially periodic state. Note that in the context of nonconserved dynamics [19, 46] one gener- 
ally selects a nonlinear term g^i directly, although this term is related to /„! through the relation 

= — (i/ni/(i0, i.e., (723 or (735. As we shall see below, in the conserved case, having the nonlinear 
term /23 describes the generic case and the role of the coefficient 62 is effectively played by the 
value of 00, which is the average value of the order parameter 0(x). 

Equation (3) can be studied in one, two or more dimensions. In one dimension with (723 the 
equation was studied by Matthews and Cox [37, 38] as an example of a system with a conserved 
order parameter; this equation is equivalent to Eq. (3) with imposed nonzero mean 0. A weakly 
localized state of the type that is of interest in the present paper is computed in [37] and discussed 
further in [38]. 
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6. Localized states in one spatial dimension 



We first consider Eq. (3) in one dimension, with a = 1 and q = 1, i.e., 

dt(P = 92 [r0+ (92 + 1)20 + 03]. (4) 

This equation is reversible in space (i.e., it is invariant under x — )■ — x). Moreover, it conserves the 
total "mass" dx, where L is the size of the system. In the following we denote the average 
value of in the system by 0o = (0) so that perturbations = — 0o necessarily satisfy (0) = 0, 
where (■ ■ ■ ) = L^^ J^{- ■ ■)dx. 

Steady states {dt = 0) are solutions of the fourth order ordinary differential equation 

= r0 + {d,, + 1)20 + 0=^ - /i, (5) 

where ji = 5F[(j)]/5(j) is an integration constant that corresponds to the chemical potential. 

Each solution of this equation corresponds to a stationary value of the underlying Helmholtz 
free energy 



F 



L 



:i + r)^ + ^-(9..0)2 + 1(9.^0)2 



dx. (6) 



We use the free energy to define the grand potential 

n = F- [ fx^dx (7) 







and will be interested in the normalized free energy density / = (F[0(x)] — F[0o])/L and in the 
density of the grand potential uj = Vl/L = F[(j){x)\/L — /i0o. We also use the LF' norm 



1 



L 



\m\ = \lj^ X (0-0o)2cia; (8) 

as a convenient measure of the amplitude of the departure of the solution from the homogeneous 
background state = 0o. 

Linearizing Eq. (4) about the steady homogeneous solution = 0o using the ansatz (50(a;, t) = 
0(a;, it) — 00 = e exp(/3t + ikx) with e ^ 1 results in the dispersion relation 

P = -e[r + {l-ef + ?,<f^. (9) 

It follows that in an infinite domain, the threshold for instability of the homogeneous state corre- 
sponds to = —300. In a domain of finite length L with periodic boundary conditions (PBC), 
the homogeneous state is linearly unstable for r < r„, where 

r„ = -(l-fc2)2_3^2^ (IQ) 
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FIG. 1: (color online) The phase diagram for the Id PFC model (3) when q = 1. The red solid lines are the 
coexistence curves between the periodic and uniform phases calculated using a two mode approximation 
[28]. The green squares show the coexistence values calculated from simulations [28]. The red circles are 
the tricritical points. The blue dashed Une shows the curve of marginal stability of the uniform state within 
linear stability theory. 

and kn = 27m/ L, n = 1,2, ... . Standard bifurcation theory with PBC shows that for L < oo 
each r„ corresponds to a bifurcation point creating a branch of periodic solutions that is uniquely 
specified by the corresponding integer n. For those integers n for which r„ > —9/2 the branch 
of periodic states bifurcates supercritically (i.e., towards smaller values of |0o|); for r„ < —9/2, 
the bifurcation is subcritical (i.e., the branch bifurcates towards larger values of \(po\)- Since each 
solution can be translated by an arbitrary amount d (mod L), each bifurcation is in fact a pitchfork 
of revolution. Although the periodic states can be computed analytically for r ^ r„, for larger 
values of |r — r„| numerical computations are necessary. In the following we use the continuation 
toolbox AUTO [47, 48] to perform these (and other) calculations. For interpreting the results it is 
helpful to think of r as a temperature-like variable which measures the undercooling of the liquid 
phase. 

Before we can discuss LS in the above model, it is helpful to refer to the phase diagram appro- 
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priate to a one-dimensional setting (Fig. 1). As shown in Ref. [28] the tricritical point is located at 
(006, r™^'') = (±a/3/38, —9/38). For r > r™^"" there exists no thermodynamic coexistence zone 
between the homogeneous and periodic states. Such a region is only present for r < r™^^ and is 
limited by the binodal lines that indicate the values of (po for which the homogeneous and periodic 
solutions at fixed r have equal chemical potential and pressure (i.e., equal grand potential). Thus 
for r < r^^"^ the transition from the homogeneous to the periodic state is of first order. The bin- 
odals can either be calculated for specific domain sizes and periods of the periodic structure or for 
an infinite domain. In the latter case the period of the periodic state is not fixed but corresponds 
to the period that minimizes the Helmholtz free energy at each (0o, r) [28]. We remark that with 
this choice of parameters the tricritical point is not the point at which the bifurcation to the pe- 
riodic state changes from supercritical to subcritical. As already mentioned, the latter occurs at 
(00, r) = (±a/3/2, —9/2), i.e., at values of r much smaller than r™^''. Further discussion of this 
point may be found in the conclusions. 

ni. RESULTS FOR THE CONSERVED SWIFT-HOHENBERG EQUATION 
A. Families of localized states 

Since Eq. (4) represents conserved gradient dynamics based on an energy functional that allows 
for a first order phase transition between the homogeneous state and a periodic patterned state, one 
may expect the existence of localized states (LS) to be the norm rather than an exception. In the 
region between the binodals, where homogeneous and periodic structures may coexist, the value 
of 00, i-C-, the amount of 'mass' in the system, determines how many peaks can form. 

As in other problems of this type we divide the LS into three classes. The first class consists of 
left-right symmetric structures with a peak in the middle. Structures of this type have an overall 
odd number of peaks and we shall refer to them as odd states, hereafter LSodd- The second class 
consists of left-right symmetric structures with a dip in the middle. Structures of this type have 
an overall even number of peaks and we refer to them as even states, hereafter LScvcn- Both types 
have even parity with respect to reflection in the center of the structure. The third class consists of 
states of no fixed parity, i.e., asymmetric states, LSasym- The asymmetric states are created from 
the symmetric states at pitchfork bifurcations and take the form of rungs on a ladder- like structure 
with interconnections between LSodd to LSeven- In view of the gradient structure of Eq. (4) the 
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asymmetric states are likewise stationary solutions of the equation. 
We now address the following questions: 

1. Do localized states exist outside the binodal region? Can they form the energetic minimum 
outside the binodal region? 

2. How does the bifurcational structure of the localized states change with changes in the 
temperature-like parameter r? How does the transition from tilted or slanted snaking to no snaking 
occur? What is the behavior of the asymmetric localized states during this process? 

Answers to these and other questions can be obtained by means of an in-depth parametric study. 
In the figures that follow we present bifurcation diagrams for localized states as a function of the 
mean order parameter value 0o for a number of values of the parameter r. All are solutions of Eq. 
(5) that satisfy periodic boundary conditions on the domain < x < L, and are characterized by 
their norm 1 |, chemical potential fi, free energy density /, and grand potential density u as 
defined in Sec. II B. 

In Fig. 2 we show the results for r = —0.9 for L = 100. Figure 2(a) shows 1 | as a function 
of 00- the classical bifurcation diagram. For these parameter values, as 0o is increased the homo- 
geneous (liquid) phase first becomes unstable to perturbations with mode number n = 16 (i.e., 16 
bumps), followed closely by bifurcations to modes with n = 15 and n = 17. All other modes 
bifurcate at yet smaller values of |0o| and are omitted. All three primary bifurcations are supercrit- 
ical. The figure also reveals that the n = 16 branch undergoes a secondary instability already at 
small amplitude; this instability creates a pair of secondary branches of spatially localized states, 
LSodd (solid line) and LScvcn (dashed line). With increasing amplitude these branches undergo 
slanted snaking as one would expect on the basis of the results for related systems with a con- 
servation law [21, 24]. The LSodd and LScvcn branches are in turn connected by ladder branches 
consisting of asymmetric states LSasym, much as in standard snaking [46]. Sample solution pro- 
files for these three types of LS are shown in Fig. 3. The snaking ceases when the LS have grown 
to fill the available domain; in the present case the LSodd and LSeven branches terminate on the 
same n = 16 branch that created them in the first place. Whether or not this is the case depends in 
general on the domain length L, as discussed further in Ref. [49]. 

The key to the bifurcation diagram shown in Fig. 2(a) is evidently the small amplitude bifurca- 
tion on the n = 16 branch. This bifurcation destabilises the n = 16 branch that would otherwise 
be stable and is a consequence of the presence of the conserved quantity 0o [37]. As L increases, 
the bifurcation moves down to smaller and smaller amplitude, so that in the limit L — )• oo the pe- 
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riodic branch is entirely unstable and the LS bifurcate directly from the homogeneous state. Since 
the LS bifurcate subcritically it follows that such states are present not only when the primary 
pattern-forming branch is supercritical but moreover are present below the onset of the primary 
instability. 

Figure 2(b) shows the corresponding plot of the free energy density / as a function of (pQ. This 
figure demonstrates that throughout much of the range of 0o the localized states have a lower free 
energy than the extended periodic states. In this range the LS are therefore energetically favored. 
Figure 2(c) shows the corresponding plot of the chemical potential /i while Fig. 2(d) shows the 
grand potential density uj. Of these Fig. 2(c) is perhaps the most interesting since it shows that the 
results of Fig. 2(a), when replotted using (0o, 1^) to characterize the solutions, in fact take the form 
of standard snaking, provided one takes the chemical potential fi as the control parameter and 0o 
as the response. In this form the bifurcation diagram gives the values of 0o that are consistent 
with a given value of the chemical potential fi (recall that 0o is related to the total particle number 
density). 

We now show how the bifurcation diagrams evolve as the temperature-like parameter r changes. 
We begin by showing the bifurcation diagrams for decreasing values of |r|. In the Appendix we 
use amplitude equations to determine the direction of branching of the localized states. Here we 
discuss the continuation results. 

The bifurcation diagram for r = —0.7 (Fig. 4(a,b)) resembles that for r = —0.9 (Fig. 2) 
although the snaking structure has moved towards smaller |0o| and is now thinner. In addition, 
it is now the second saddle-node on the LSodd branch that lies farthest to the left, and not the 
first. For r = —0.6 (Fig. 4(c,d)), the branches of localized states still form a tilted snakes-and- 
ladders structure, but the saddle-nodes on the LSodd and LScvcn branches are now absent, i.e., 
both solution branches now grow monotonically. The resulting diagram has been called "smooth 
snaking" [50]. However, despite the absence of the saddle-nodes on the LSodd and LSeven branches 
the interconnecting ladder states consisting of asymmetric states still remain. This continues to be 
the case when r = —0.5 (Fig. 4(e,f)) although the structure has moved to yet larger 0o and the 
snake has become even thinner. Finally, for r = —0.4 (Fig. 4(g,h)), the snake is nearly dead, and 
only tiny wiggles remain. The bifurcation of the localized states from the n = 16 periodic state is 
now supercritical (see Appendix) but the LS branches continue to terminate on the same branch at 
larger amplitude, and do so via a single saddle-node at the right (Fig. 4(h)). Sample profiles along 
the resulting LSodd branches are shown for several values of r in Fig. 5. We note that a change of r 
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(C) ^0 t (d) 

FIG. 2: (color online) Characteristics of steady state (localized) solutions of the one-dimensional conserved 
Swift-Hohenberg equation (4) as a function of the mean order parameter (po for a fixed domain size of 
L = 100 and r = —0.9. The various solution profiles are characterized by their (a) norm ||5</>||, (b) 
free energy density /, (c) chemical potential fi, and (d) grand potential density u. The thick green dotted 
line, labeled 'flat', corresponds to the homogeneous solution (j){x) = (pQ. Periodic solutions with n = 16 
peaks are shown as a thin green dashed line, whereas the nearby thin blue and black dotted lines represent 
the n = 15 and n = 17 solutions, respectively. The thick solid black and dashed red lines that bifurcate 
from the n = 16 periodic solution represent symmetric localized states with a maximum (LSodd) and a 
minimum (LSeven) at their center, respectively. Both terminate on the n = 16 solution. The 14 blue solid 
lines that connect the LSodd and LSeven branches of symmetric localized states correspond to asymmetric 
localized states (LSasym)- Together these three sets of branches of localized states form a tilted snakes-and- 
ladders structure. Typical order parameter profiles along the three LS branches are shown in Fig. 3, and 
correspond to locations indicated in panels (a) and (c) by filled black squares (LSodd)> red circles (LSeven). 
blue triangles (LSasym) and green diamond (periodic solution with n = 16). 
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FIG. 3: (Color online) A selection of steady state profiles (j){x) — (j)Q for r = —0.9 and values of 4>q in 
the range —0.65 < (j)o < —0.3 (the number in each panel indicates the corresponding value of ^q). Going 
from top left to bottom right we first show nine LSodd solutions, i.e., symmetric localized states with an 
odd number of maxima (in black), then eight LSoven solutions, i.e., symmetric localized states with an even 
number of maxima (in red), followed by six LSasym solutions, i.e., asymmetric localized states (in green). 
The final plot is the n = 16 periodic solution for (/iq = —0.3 (in blue). The solutions on the symmetric 
branches correspond to locations indicated in Fig. 2(a) and are shown in order, starting from the bifurcation 
point that creates them and continuing to larger norm 1 1(5</>| |. The color coding corresponds to that used in 
Figs. 2(a,c): LSodd (filled black squares), LSeven (red circles), LSasym (blue triangles) and periodic (green 
diamond). 

has a profound effect on the transition region between the homogeneous background state and the 
periodic state: with decreasing |r| the LS become wider and the localized periodic structure looks 
more and more like a wave packet with a smooth sinusoidal modulation of the peak amplitude. 

As the 'temperature' r decreases even further, the bifurcation diagrams remain similar to those 
displayed in Fig. 2 until just before r = —1.5, where substantial changes take place and the 
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FIG. 4: (color online) The norm (left) and chemical potential (right) of the homogeneous, periodic and 
localized steady state solutions as a function of the mean order parameter 0o, for a fixed domain size of 
L = 100 and several values of r > —0.9. In (a,b) r = —0.7, (c,d) r = —0.6, (e,f) r = —0.5, and (g,h) 
r = —0.4. The line styles are as in Fig. 2. Typical order parameter profiles along the branches of symmetric 
localized states with an odd number of maxima (black lines) are shown in Fig. 5, and correspond to locations 
indicated in the panels by filled black squares. 13 
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FIG. 5: (Color online) A selection of steady state profiles (j){x) — along the LSodd branches shown in 
Fig. 4 for r = —0.7 (1st row), r = —0.6 (2nd row), r = —0.5 (3rd row), r = —0.4 (4th row), and 
r = —0.375 (last row), for various values of 0o> as indicated in the top left corner of each panel. The 
solutions in each row are shown in order, starting from near the bifurcation point that creates them and 
continuing to larger norm \\S(p\\. The locations of the profiles are indicated in Fig. 4 by the filled black 
squares. The bifurcation diagram for r = —0.375 (not shown in Fig. 4) is qualitatively the same as that for 
r = -0.4. 

complexity of the bifurcation diagram grows dramatically. This is a consequence of the appearance 
of other types of localized states that we do not discuss here. Likewise, we omit here all bound 
states of the LS described above. These are normally found on an infinite stack of isolas that are 
also present in the snaking region [51, 52]. 



B. Tracking the snake 

In Fig. 6 we show, for — 1.5 < r < —0.4, the result of tracking all the saddle-node bifurcations 
visible in the previous bifurcation diagrams in the (0o, r) plane, while Fig. 7 shows an enlarge- 
ment of the region —1 < r < —0.35 together with the result of tracking the tertiary pitchfork 
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FIG. 6: (color online) Loci of saddle-node bifurcations on the branches of symmetric localized states in 
the (005 plane for r > —1.5. Saddle-nodes annihilate pairwise as r increases (solid black lines: LSodd; 
dashed red lines: LSevcn)- 

bifurcations to the asymmetric states. 

Figure 6 shows that the saddle-nodes annihilate pairwise in cusps as r increases. The annihila- 
tions occur first for smaller 0o and later for larger (f)Q, and occur alternately on LSodd and LSgven- 
Above the locus of the cusps the snaking is smooth, although as shown in Fig. 7 the tertiary ladder 
states remain. The thick green curve in Fig. 7 represents the locus of the secondary bifurcation 
from the n = 16 periodic state to LS and shows that on either side the bifurcation to LS is subcrit- 
ical for sufficiently negative r but becomes supercritical at larger r, cf. Fig. 4(h) and Appendix. 

One may distinguish six intervals in r with different types of behavior. These depend on the 
system size, but the results in Fig. 7 for L = 100 are representative: 

(i) For r > —0.33: no LS exist, the only nontrivial states are periodic solutions. 

(ii) For —0.33 > r > —0.39: branches of even and odd symmetric LS are present, and appear 
and disappear via supercritical secondary bifurcations from the branch of periodic solutions. 
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FIG. 7: (color online) Loci of the saddle-node bifurcations on the branches of symmetric localized states 
in the {(pQ, r) plane (solid black lines: LSodd; dashed red lines: LSevcn) together with the bifurcations to 
asymmetric localized states (dotted blue lines) for r > —1.5. The thick green curve represents the locus of 
secondary bifurcation from the periodic n = 16 state to localized states and shows that on either side the 
bifurcation to LS is subcritical for sufficiently negative r but becomes supercritical at larger r. 

With decreasing r, more and more branches of asymmetric LS emerge from these two sec- 
ondary bifurcation points. 

(iii) For —0.39 > r > —0.41: both branches of symmetric LS emerge subcritically at large (po 
and supercritically at small 0o- 

(iv) For —0.41 > r > —0.56: both branches of symmetric LS emerge subcritically at either end. 
Further branches of asymmetric LS emerge with decreasing r either from the two secondary 
bifurcation points or from the saddle-node bifurcations on the branches of symmetric LS, 
but the symmetric LS still do not exhibit snaking, i.e., no additional folds are present on the 
branches of symmetric LS. 
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(v) For —0.56 > r > —0.64 (highlighted by the grey shading in Fig. 7): pairs of saddle-nodes 
appear successively in cusps as r decreases, starting at larger 0o- Thereafter saddle-nodes 
appear alternately on branches of even and odd symmetric LS. The appearance of the cusps 
is therefore associated with the transition from smooth snaking to slanted snaking. 

(vi) For —0.64 > r: The slanted snake is fully developed. Only one further pair of saddle-node 
bifurcations appears in the parameter region shown in Fig. 7. With decreasing r the snaking 
becomes stronger; each line in Figs. 6 and 7 that represents a saddle-node bifurcation crosses 
more and more other such lines, i.e., more and more different states are possible at the same 
values of 0o- Furthermore, the subcritical regions (outside the green curve in Fig. 7) become 
larger. 

C. Relation to binodal lines 

From the condensed matter point of view, where the cSH/PFC equation represents a model for 
the liquid (homogeneous) and solid (periodic) phases, one is particularly interested in results in 
the thermodynamic limit L — )• oo. As mentioned above in the context of the phase diagram in 
Fig. 1, the binodal lines correspond to values of (0o, r) at which the homogeneous state and the 
minimum energy periodic state coexist in the thermodynamic limit. These are defined as pairs of 
points at which the homogeneous state and the periodic state have the same 'temperature' (i.e., 
same r value), the same chemical potential fi and the same pressure p = —cu, and are displayed as 
the blue dash-dot lines in Fig. 8. For a given value of r, these two lines give the values of 0o of the 
coexisting homogeneous (lower 0o) and periodic (higher 0o) states. Note that when plotted with 
the resolution of Fig. 8, the binodals are indistinguishable from the coexistence lines between the 
finite size L = 100, n = 16 periodic solution and the homogeneous state. Figure 8 also displays 
the line (green solid line) at which the L = 100 localized states bifurcate from the n = 16 branch 
of periodic solutions. For r < — 0.4 the LS bifurcations are subcritical implying that the localized 
states are present outside the green solid line. Figure 8 shows the loci of the outermost saddle-node 
bifurcations on the branches of symmetric LS that result (dashed lines to the left and to the right 
of the green solid line for r < —0.4). The most striking aspect of Fig. 8 is that for r < — 1 these 
Unes actually cross and exit the region between the two binodal lines, indicating that in the PFC 
model one can find stable LS outside of the binodal. Although these are not the lowest free energy 
states (we have checked this for r > —1.5), this remarkable fact points towards the possibility 
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FIG. 8: (color online) Loci in the {(pQ, r) plane computed for L = 100 of (i) the primary bifurcation from the 
homogeneous state to the n = 16 periodic branch (dotted orange line); (ii) the bifurcation of the localized 
states from the n = 16 periodic branch (solid green line); (iii) the outermost saddle-node bifurcations on the 
branches of symmetric localized states (long-dashed red and short-dashed black lines); and (iv) the binodals 
between the periodic and homogeneous states (dot-dashed blue lines). The latter coincide with the binodals 
for the n = 16 periodic state and the homogeneous state to within the resolution of the figure. Figures 6 and 
7 show how the loci of the outermost saddle-node bifurcations of localized states fit into the overall picture. 

of metastable nanocrystals existing outside of the binodal. We must mention, however, that these 
structures have been found in a finite size system with L = 100; we have not investigated their 
properties for larger system sizes L. 

D, Stability 

Although we have not computed the stability properties of the different localized states with 
respect to infinitesimal perturbations, we can use the principles of bifurcation theory to deduce the 

18 

Prenrint- contact: u.thiele@lhoro.ac.uk - www.uwethiele.de - January 21 . 201 3 



likely stability properties. In the reference bifurcation diagram in Fig. 2(a) the homogeneous state 
(liquid) is stable for large negative (po, and for a system with L = 100 loses stability to the n = 16 
mode as 0o increases. Since the bifurcation is supercritical, the n = 16 periodic state is initially 
stable. The LS that bifurcate from it subcritically will both be unstable. The single bump state 
is likely to be once unstable and it therefore acquires stability at the first saddle-node on the left, 
remaining stable until the first crosslink where a second (phase) eigenvalue becomes unstable. The 
first (amplitude) eigenvalue becomes unstable at the saddle-node on the right so that the portion of 
LSodd with negative slope is twice unstable. These instabilities are then undone so that the LSodd 
above the second crosslink are again stable. The transition to smooth snaking that occurs with 
decreasing r eliminates instabilities associated with the first eigenvalue but not the second. 

The LScvcn states are initially twice unstable but both eigenvalues stabilize near the first saddle- 
node on the left, so that LScvcn is stable on the part of the branch with positive slope, below 
the second crosslink etc. Thus the connecting LSasym are unstable. Once again, the transition to 
smooth snaking that occurs with decreasing r eliminates instabilities associated with the amplitude 
eigenvalue but not those associated with the phase eigenvalue. One can check that in this case the 
connecting LSasym remain unstable with the LSodd stable between the second and third crosslinks 
(Fig. 4(c)), and unstable between the third and fourth crosslinks. Likewise the LScvcn are stable 
above the left saddle-node and below the second crosslink; they are unstable between the second 
and third crosslinks and acquire stability at the third crosslink etc (Fig. 4(c)). Note that as a result 
of these stability assignments there is at least one LS state that is stable at all 0o values between 
the leftmost and rightmost saddle-nodes. 

IV. LOCALIZED STATES IN TWO AND THREE DIMENSIONS 
A. Numerical algorithm 

To perform direct numerical simulations (DNS) of the conserved Swift- Hohenberg equation in 
higher dimensions we use a recently proposed algorithm [53] that has been proved to be uncondi- 
tionally energy-stable. As a consequence, the algorithm produces free-energy-decreasing discrete 
solutions, irrespective of the time step and the mesh size, thereby respecting the thermodynamics 
of the model even for coarse discretizations. For the spatial discretization we employ Isogeomet- 
ric Analysis [54], which is a generalization of the Finite Element Method. The key idea behind 
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't'o 

FIG. 9: (color online) Phase diagram for the conserved Swift-Hohenberg equation, Eq. (3), in two dimen- 
sions when q = I. The red solid lines show the various coexistence curves, the blue dotted line shows the 
limit of linear stability and the grey striped areas show phase coexistence regions. The small panels show 
steady solutions at r = —0.9 and selected mean order parameter values 0o indicated by the black squares. 
From top to bottom left (j)o = —0.15, —0.2 and —0.25, while from top to bottom right (j)o = —0.3, —0.35 
and —0.725. Additional solutions that detail the transition occurring between the last two panels are shown 
in Fig. 10. The domain size is 100 x 100. 

Isogeometric Analysis is the use of Non-Uniform Rational B-Splines (NURBS) instead of the 
standard piecewise polynomials used in the Finite Element Method. With NURBS Isogeometric 
Analysis gains several advantages over the Finite Element Method. In the context of the con- 
served Swift-Hohenberg equation, the most relevant one is that the Isogeometric Analysis permits 
the generation of arbitrarily smooth basis functions that lead to a straightforward discretization of 
the higher-order partial derivatives of the conserved Swift-Hohenberg equation [55]. For the time 
discretization we use an algorithm especially designed for the conserved Swift-Hohenberg equa- 
tion. It may be thought of as a second-order perturbation of the trapezoidal rule which achieves 
unconditional stability, in contrast with the trapezoidal scheme. All details about the numerical 
algorithms may be found in [53]. 
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B. Two dimensions 

As in one spatial dimension, the phase diagram for two-dimensional structures helps us to 
identify suitable parameter values where LS are likely to occur. The phase diagram in Fig. 9, 
determined numerically [28], shows three distinct phases, labeled 'bumps', 'stripes' and 'holes'. 
In view of the fact that = 0o + S(f){x), bumps and holes correspond to perturbations (5</)(x) 
with opposite signs; both have hexagonal coordination. In addition, the phase diagram reveals four 
regions of thermodynamic coexistence (hatched in Fig. 9), between bumps and the uniform state, 
between bumps and stripes, between stripes and holes, and between holes and the uniform state, 
respectively. 

Examples of results obtained by DNS of Eq. (3) starting from random initial conditions are 
displayed in the side panels of Fig. 9. These six panels show results for fixed r = —0.9. For small 
values of 1 00 1 the system forms a labyrinthine lamellar-like stripe state; the stripes pinch off locally 
into bumps as |0o| is increased, leading to the formation of inclusions of bumps in a background 
stripe state. The pinching tends to occur first at the ends of a stripe and then proceeds gradually 
inwards. In other cases, free ends are created by the splitting of a stripe into two in a region of 
high curvature. The formation of bumps tends to take place at grain boundaries, and once bump 
formation starts, it tends to spread outward from the initial site. For 0o = —0.2 the areas covered 
by stripes and bumps are comparable and for larger values of |0o| the bump state dominates. By 
00 = —0.3 the stripes are almost entirely gone and the state takes the form of a crystalline solid 
with hexagonal coordination but having numerous defects. As |0o| increases further, vacancies 
appear in the solid matrix and for large enough |0o| the solid "melts" into individual bumps or 
smaller clusters, as further detailed in Fig. 10. 

Figure 10 shows further results from a scan through decreasing values of 0o at r = —0.9. 
We focus on the relatively small range 0o = —0.45 to 0o = —0.675, where localized states 
occur. These reveal a gradual transition from a densely packed solid-like structure to states with a 
progressively increasing domain area that is free of bumps, i.e., containing the homogeneous state. 
The bumps percolate through the domain until approximately 0o = —0.575. For smaller values 
of 00 the order parameter profiles resembles a suspension of solid fragments in a liquid phase; the 
solid is no longer connected. As 0o decreases further the characteristic size of the solid fragments 
decreases, as the solid fraction falls. 
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FIG. 10: (color online) Steady-state solutions of the conserved Swift-Hohenberg equation, Eq. (3), in two 
dimensions for r = —0.9 and different values of in the range —0.675 < 00 < —0.45, where localized 
states occur. The corresponding value of is indicated below each panel. The domain size is 100 x 100. 



C. Three dimensions 

In three dimensions, Eq. (3) exhibits a large number of steady state spatially periodic structures. 
These include those with the symmetries of the simple cubic lattice, the face-centered cubic lattice 
and the body-centered cubic lattice [56]. Although we do not calculate the phase diagram for 
the three-dimensional (3d) system, numerical simulations in 3d reveal that a lamellar (parallel 
'sheets') state is energetically preferred for small 0o- Slices through these structures resemble the 
stripes observed in 2d. As |0o| increases, the lamellae pinch off, much as in two dimensions, and 
progressively generate a 3d disordered array of bumps (Fig. 11). This solid-like state is far from 
being a perfect crystal, however, and with increased |0o| develops vacancies which eventually lead 
to its dissolution, just as in two dimensions. 

In Fig. 12 we superpose the 2d and 3d DNS results for r = —0.9 on the Id bifurcation diagrams 
for r = —0.9 and system size L = 100 (Fig. 2). The 2d order parameter profiles are calculated 
on a domain with area = 100^ and the 3d results have a system volume = 100^. In Fig. 
12(a) we display the norm in (b) the chemical potential /i, in (c) the mean free energy 

(F — Fq)/L'^, and in (d) the mean grand potential to = F/L'^ — (f)Qjj,. In each plot the connected 
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FIG. 11: (color online) Steady-state localized solutions of Eq. (3) in three dimensions for r = 
—0.9 and different mean order parameter values (pQ-. from top left to bottom right 0o = 
0.025, 0.125, 0.225, 0.325, 0.425, 0.525, 0.625, 0.725 and 0.750. The domain size is 100 x 100 x 100. 

(violet) squares correspond to results from 2d calculations such as those displayed in Figs. 9 and 
10 and the connected green triangles correspond to 3d results such as those displayed in Fig. 11. 
It is from examining the figures in panels (c) and (d) that one can most easily discern the reason 
for the main differences between the Id results and the 2d and 3d results: we see (particularly in 
the 2d results) that the chemical potential ^ and the pressure p = -co have regions where these 
measures are roughly flat as a function of 0o and regions where they increase as a function of 0o- 

23 

Prenrint- contact: u.thiele@lboro.ac.uk - www.uwethiele.de - January 21 . 201 3 



We also see that these increases in some places relate to features in the Id results but in other 
places they do not have any relation to what one sees from the Id results. This is because in 2d and 
3d the system displays phases that are not seen in Id (cf. Figs. 1 and 9). To understand the origin 
of these roughly flat portions, we recall that in the thermodynamic limit L oo two states are 
said to be at coexistence if the 'temperature' r, the chemical potential fi and the pressure p = —uj 
is the same for both states. These quantities do not change in value as one takes the system across 
the coexistence region by increasing the average density in the system (or equivalently 0o)- This 
is because the additional surface excess free energy terms (surface tension terms that are present 
because both coexisting phases are in the system) do not contribute in the thermodynamic hmit. 
This is in turn a consequence of the fact that these surface terms scale as L^'^~^\ whereas the 
bulk volume terms scale as L^, where d is the dimensionality of the system. Thus, regions where 
these measures are approximately flat are in a coexistence region between two phases - this can 
be confirmed for the 2d results by comparing the ranges of (pQ, where the results for fi and cu are 
approximately flat, with the coexistence regions in Fig. 9. The observation that the 2d curves in 
Fig. 12 are not completely flat indicates that the interfacial (surface tension) terms between the 
different phases in the LS state do contribute to the free energy, and is thus a finite size effect. 
Note that it might be possible to distribute the LS in such a way that they percolate throughout the 
whole system so that the contribution from the interfaces scales as L'^. If this is the case, the above 
argument does not apply. 



V. DISCUSSION AND CONCLUSIONS 



The conserved Swift-Hohenberg equation is perhaps the simplest example of a pattern-forming 
system with a conserved quantity. Models of this type arise when modelling a number of different 
systems, with the PFC model being one particular example. Other examples include binary fluid 
convection between thermally insulating boundaries [26], where this equation was first derived, 
convection in an imposed magnetic field (where the conserved quantity is the magnetic flux [20, 
21]) and two-dimensional convection in a rotating layer with stress-free boundaries (where the 
conserved quantity is the zonal velocity [20, 22]). Models of a vibrating layer of granular material 
are also of this type (here the conserved quantity is the total mass [50]). It is perhaps remarkable 
that all these distinct systems behave very similarly. In particular, they all share the following 
features: (i) strongly subcritical bifurcations forming localized structures so that the resulting 
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(C) ^0 ^0 (d) 

FIG. 12: (color online) Characteristics of steady-state localized solutions of the conserved Swift-Hohenberg 
equation for r = —0.9, as a function of the mean order parameter 0o on fixed one, two and three- 
dimensional domains of size L"', d=l,2, 3, and L = 100. The various solution profiles are characterized 
by their (a) norm 1 1 | , (b) chemical potential fi, (c) mean free energy {F — Fq)/ L'^, and (d) mean grand 
potential u = F/L'^ — (pQjjL. The connected (violet) square and (green) triangle symbols correspond to 2d 
calculations (for sample profiles see Figs. 9 and 10) and 3d calculations (for sample profiles see Fig. 11). 
For comparison we also show the Id results from Fig. 2 for the periodic (i.e., stripe) state with n = 16 
bumps, green dashed) and the various Id localized states. The thick green dotted line corresponds to the 
homogeneous solution i;^(x) = ^q. 

LS are present outside of the bistability region between the homogeneous and periodic states; 
(ii) presence of LS even when the periodic branch is supercritical; (iii) organization of LS into 
slanted snaking; and (iv) the transition from slanted snaking to smooth snaking whereby the LS 
grow smoothly without specific bump-forming events (referred to as 'nucleation' events in the 
pattern formation literature [84].) These properties of the system can all be traced to the fact that 
the conserved quantity is necessarily redistributed when an instability takes place or a localized 
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FIG. 13: (color online) The norm \ \6(j)\ \ for the homogeneous, periodic and localized steady state solu- 
tions of the conserved Swift-Hohenberg equation (3) as a function of the chemical potential fi, for a fixed 
domain size of L = 100 and various r. The horizontal thick dot-dashed green line corresponds to the ho- 
mogeneous solution (j){x) = (pQ. Periodic solutions with n = 16 peaks are labeled by the corresponding 
r values, whereas the branches that bifurcate from the periodic solution branches represent the two types 
of symmetric localized states (LSodd and LSeven): f = —0.7 (dotted red lines), r = —0.6 (dashed black 
lines), r = —0.5 (solid blue lines), r = —0.4 (dotted green lines), r = —0.3 (solid black line) and r = —0.2 
(dashed red line). The dashed vertical lines indicate the corresponding coexistence chemical potential values 
for r = —0.7, —0.6 and —0.3, respectively. 

structure forms. This fact makes it harder for additional LS to form and as a result the system 
has to be driven harder for this to occur, leading to slanted snaking. Bistability is no longer 
required since the localized structures are no longer viewed as inclusions of a periodic state within 
a homogeneous background or vice versa. These considerations also explain why the bifurcation 
diagrams in these systems are sensitive to the domain size, and it may be instructive, although 
difficult, to repeat some of our calculations for larger domain sizes. 

Of particular significance is the observation that the slanted snaking found when the LS norm 
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\\S(f)\\ is plotted as a function of the mean order parameter 0o for fixed r (Figs. 2 and 4(a,b)) is 
'straightened' out when the same solution branches are displayed as a function of yU (see Fig. 13). 
The snaking is then vertically aligned and centred around the coexistence chemical potential values 
(referred to as Maxwell points in Ref. [46]) as known from the non-conserved Swift-Hohenberg 
equation [46]. In this representation the periodic states typically bifurcate subcritically from the 
homogeneous state (Fig. 13), in contrast to the supercritical transitions found when 0o is used as 
the bifurcation parameter (Figs. 2 and 4). It is also worth noting that the thermodynamic tricritical 
point at (006, r^^"^) = (±a/3/38, —9/38) discussed in Sec. II B (see also the last paragraph of the 
Appendix) corresponds to the transition from a subcritical to a supercritical primary bifurcation in 
the representation of Fig. 13, i.e., when yU is used as the bifurcation parameter. This transition takes 
place between the solid black line for r = —0.3 and the dashed red line for r = —0.2 in Fig. 13 
and implies that the linear stability properties of identical periodic solutions differ, depending on 
whether the permitted perturbations preserve /i (and hence permit 0o to vary) or vice versa, and 
this is so for the localized states as well. In particular, when fi is used as the control parameter, 
the differences identified in Sec. IIIB between cases (i)-(iv) where the LS branches do not snake, 
and cases (v)-(vi) where they do snake disappear. In Fig. 13 snaking over n is clearly visible for 
r = —0.7, —0.6 and —0.5. There is no snaking for r = —0.4, most likely due to finite size effects, 
or for larger r. 

The thermodynamic reason that the snaking becomes straightened when displayed as a function 
of the chemical potential /i is related to the issues discussed at the end of Sec. IV. For a system with 
(00, r) chosen so that it is in the coexistence region (e.g., between the homogeneous and the bump 
states) and size L large enough to be considered to be in the thermodynamic limit L — )■ oo, the 
chemical potential does not vary when 0o is changed and neither does the grand potential cu = —p. 
This is because the free energy contribution from interfaces between the coexisting phases (the 
interfacial tension) scales with the system size ~ L^'^^^\ and hence is negligible compared to the 
bulk contributions which scale ~ L"'. As 0o is varied so as to traverse the coexistence region, new 
bumps are added or removed from the bump state (LS). For a (thermodynamically) large system, 
the resulting changes in the free energy are negligible, although this is no longer true of a finite 
size system. Furthermore, the interfacial free energy contribution between the two phases varies 
depending on the size of the bumps right at the interface and the size of these bumps depends on 
the value of 0o- The difference between the maximal interfacial energy and the minimal value 
determines the width of the snake. All of these contributions have leading order terms that scale ~ 
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L . In a related situation, the variation in the chemical potential as the density is varied in a finite- 
size system containing a fluid exhibiting gas-liquid coexistence is discussed in Refs. [57-59]. It is 
instructive to compare, for example. Fig. 3 in Ref. [59] and Fig. 12(c) of the present paper. 

We mention that the steady states of the nonconserved Swift-Hohenberg equation studied, e.g., 
in Refs. [46, 49], correspond to solutions of Eq. (5) with the nonlinearity (723 but yU = (see 
Sec. II A). These always show vertically aligned snaking when the LS norm \ \6(f)\ \ is plotted as a 
function of r (adapting 0o)- However, when = no localized states are present with the pure 
cubic nonlinearity employed in Eq. (5). To find such states fi must be fixed at a value sufficiently 
far from zero; the resulting LS then exhibit vertically aligned snaking when plotted as a function 
ofr [46]. 

The domain size we have used, L = 100, is moderately large. It contains of the order of 16 
wavelengths of the primary structure-forming instability. Because the equation is simpler than the 
hydrodynamic equations for which similar behavior was observed, the results we have been able 
to obtain are substantially more complete, even in one dimension, than was possible elsewhere. In 
particular, we have been able to compute the rungs of asymmetric LS and to study their behavior 
as the system transitions from slanted snaking to smooth snaking as r increases towards zero. The 
extension of our results to two and three spatial dimensions is necessarily incomplete, although 
the transition to clusters of bumps followed by isolated bumps as the total "mass" decreases is 
not surprising. However, the transition from a connected structure to a disconnected one (the 
"percolation" threshold) in two and three dimensions deserves a much more detailed study than 
we have been able to provide. 

In this connection we mention two experimental systems exhibiting a transition from a solid- 
like phase to a gas-like phase of individual spots. This is the gaseous discharge system studied by 
H.-G. Purwins and colleagues [17, 60] and the liquid crystal light valve experiment of S. Residori 
and colleagues [2]. In both these systems a crystal-like structure of spots with hexagonal coor- 
dination was observed to melt into a 'gas' of individual spots as a parameter was varied. This 
two-dimensional process leads to states resembling those found here in Figs. 9 and 10 although 
the stripe-like structures were typically absent. In these two systems the spots in the 'gas-phase' 
are mobile unlike in the cSH equation indicating absence of variational structure. However, both 
systems are globally coupled, by the imposed potential difference in the discharge system and the 
feedback loop in the liquid crystal light valve experiment, raising the possibility that the global 
coupling in these systems plays a similar role to the role played by the conserved order parameter 
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in the cSH equation. 

We should also mention that some of the localized states observed in Figs. 10 and 11 raise 
some concern about the validity of the PFC as a model for solidification and freezing - we refer 
in particular to the order parameter profiles displayed at the bottom right of these figures. These 
show that in both 2d and 3d the PFC predicts the existence of steady states with isolated single 
bumps. Recall that in the standard interpretation of the PFC, the bumps correspond to frozen 
particles whilst the homogeneous state corresponds to the uniform liquid state. Such profiles 
could perhaps be a signature of the dynamical heterogeneity that is a feature of glassy systems, 
but there are problems with this interpretation - the glass transition is a collective phenomenon 
- single particles do not freeze on their own whilst the remainder of the particles around them 
remain fluid! We refer readers interested in the issue of the precise interpretation of the order 
parameter in the PFC to the discussion in the final section of Ref. [61]. We should also point out 
that although these structures correspond to local minima of the free energy (i.e., they are stable) 
they do not correspond to the global minimum. These states occur for state points where the global 
free energy minimum corresponds to the uniform homogeneous state. 

All the results presented here have been obtained for the generic conserved Swift-Hohenberg 
equation, Eq. (1), with the energy functional in Eq. (2). We believe that our main results provide a 
qualitative description of a number of related models in material science that are of a similar struc- 
ture and describe systems that may show transitions between homogeneous and patterned states 
characterized by a finite structure length. In particular, we refer to systems that can be described 
by conserved gradient dynamics based on an underlying energy functional that features a local 
double-well contribution, a destabilizing squared gradient term and a stabilizing squared Lapla- 
cian term. The latter two terms may themselves result from a gradient expansion of an integral 
describing nonlocal interactions such as that required to reduce DDFT models to the simpler PFC 
model [27, 28]. 

Other systems, where the present results may shed some light, include diblock copolymers and 
epitaxial layers. The time evolution equation for diblock copolymers [62-66] is of fourth order 
like the nonconserved SH equation but contains a global coupling term that is related to mass 
conservation. The equation emerges from a nonlocal term in the energy functional [67, 68]. The 
global coupling results in an evolution towards a state with a given mean value for the density 
order parameter (pQ if the initial value is different from 0o or in a conservation of mass as the 
system evolves if the initial value coincides with the imposed 0o- Although this differs from 
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the formulation using a conserved Swift-Hohenberg equation, the steady versions of the diblock- 
copolymer equation and of the conserved Swift-Hohenberg equation are rather similar: they only 
differ in the position of the nonlinearity. Up to now no systematic study of localized states exists 
for the diblock-copolymer equation, although Ref. [66] discusses their existence and gives some 
numerical examples for a profile with a single bump in rather small systems (see their Fig. 5). Since 
in the diblock-copolymer system the order parameter is a conserved quantity, we would expect the 
snaking of localized states for fixed 00 to be slanted similar to our Fig. 2, instead of being vertical, 
corresponding to a standard snake, where all the saddle-node bifurcations are vertically aligned, 
as sketched in Fig. 6 of [66]. 

Finally, we briefly mention a group of model equations that are derived to describe the evolution 
of the surface profile of epitaxially strained solid films including, e.g., the self-organization of 
quantum dots [69]. The various evolution equations that have been employed account for the 
elasticity (linear and nonlinear isotropic elasticity, as well as misfit strain) of the epitaxial layer and 
the (isotropic or anisotropic) wetting interaction between the surface layer and the solid beneath. 
The evolution equations we wish to highlight are of sixth order [70-73] much like the conserved 
Swift-Hohenberg equation investigated here. Other models, however, are of fourth order only [74, 
75] or contain fully nonlocal terms (resulting in integro-differential equations) [76, 77]. However, 
even the sixth order models often contain additional nonlinear terms in the derivatives (see, for 
instance, Eq. (5) of Ref. [69]). Localized state solutions of these equations have to our knowledge 
not yet been studied systematically, although some have been obtained numerically (Fig. 5 of 
[69]). Future research should investigate how the characteristics of the localized states analysed 
here for the conserved Swift-Hohenberg equation differ from those in specific applied systems 
such as diblock copolymers or epitaxial layers. 

Appendix 

In this Appendix we determine the direction of branching of the localized states when they 
bifurcate from the branch of periodic states. When the domain is large this bifurcation occurs 
when the amplitude of the periodic states is small and hence is accessible to weakly nonlinear 
theory. 

We begin with Eq. (4) which may be written 

0, = adl[{r + g^)0 + <f + 2q^dl<\> + 5^0]. (11) 
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This equation has the homogeneous solution cj) = (pQ. We let = 0o + obtaining 

= adllir + q^ + 30^)^ + 30o^' + + 2q^dlij + d^ij]. (12) 

Linearizing and looking for solutions of the form ip oc exp(/3t + ikx) we obtain the dispersion 
relation 

/3 = _«A;2[r+(g2_P)2 + 302] (^3) 
The condition /3 = gives the critical wavenumbers 



ki = q'±^-r-3<Pi (14) 

Thus when r = rc = — 30o the neutral curve has maxima at both k = and kc = q. When 
r = — 30Q — e^z/, where u = 0{\) and e is a small parameter that defines how far r is from then 
a band of wavenumbers near k = q grows slowly with growth rate (3 = C(e^), while wavenumbers 
near k = decay at the same rate. There is therefore time for these two disparate wavenumbers to 
interact and it is this interaction that determines the direction of branching. 

These considerations suggest that we perform a two-scale analysis with a short scale x = 
0{q^^) and a long scale X = ex, so that dx dx + edx etc. We also write 

^ = eA{X, t)e'«" + e^B{X, t) + e^C{X, t)e^"i'' + c.c. + C(e=^), (15) 

where the amplitudes A and C are complex and B is real. Substituting, we obtain 

At = -e^aq^-uA - Aq^Axx + 60o^5 + Q(j)oCA* + 3\A\'^A) + 0{e^), (16) 

Bt = e^a{q*Bxx + QM^lxx) + 0{e'), (17) 

and 

Ct = -Aaq\9q^C + 300^^) + 0(e). (18) 

The last equation implies that the mode C decays on an 0{1) time scale to its asymptotic value, 
C = —(poA^/^q^ + C(e). The resulting equations may be written in the form 

A = uA + AAxx - iAOx - 3 (l - \A\^A + 0(e), (19) 

9t = exx + m\ + 0{t) (20) 

using the substitution B = 9x, and integrating Eq. (17) once to obtain an equation for 9. In writing 
these equations we have absorbed q into the length scale X and e'^aq^ into the time scale t, and 

31 

Prenrint- contact: u.thiele@lboro.ac.uk - www.uwethiele.de - January 21 . 201 3 



introduced the parameter ^ = Qcly^/q^ < 0. The resulting equation are equivalent to the equations 
studied by Matthews and Cox [37]. 

Equations (19)-(20) provide a complete description of the small amplitude behavior of Eq. (4). 
The equations inherit a gradient structure from Eq. (4), 

where 

F[A, A\ 9] = jJy-M? + Mx? + \el + i\A\^dx + ^ (l - fj) 1^1'} d^- (22) 
We may write this energy in the form 

F[A,A\e] = jJ^-u\A\' + A\Ax? + \ (^x + ^l^lT + 3 (l - ^) 1^1^ | rfX, (23) 

implying that the free energy F[A^A* ^6] is not bounded from below once > 54/19 (equiv- 
alently cfl > 3g^/38). This is a reflection of the presence of subcritical branches. Indeed, the 
steady states of Eqs. (19)-(20) correspond to critical points of this energy and satisfy the nonlocal 
equation 

4Axx + - {\A\'))A - 3(^1 - \A\'A = 0, (24) 

where ((■■■)) = {■ ■ ■) dx and L is the domain length. This equation demonstrates (i) that 

the primary bifurcation is subcritical when > 54/19 in agreement with Eq. (23), and (ii) that 
as I A| increases the value of u has to be raised in order to maintain the same value of the effective 
bifurcation parameter z/efr = u — {\A\'^). This is the basic reason behind the slanted structure in 
Fig. 2(a). Equations of this type arise in numerous applications [22, 78-80] and their properties 
have been studied in several papers [79, 81-83]. We mention, in particular, that unmodulated 
wavetrains bifurcate supercritically when < 54, a condition that differs from the corresponding 
condition < 54/19 for spatially modulated wavetrains. 

Equations (19)-(20) possess the solution {A,d) = {Ao,0), where \Aq\'^ = u/[3{l - ^V54)], 
corresponding to a spatially uniform wavetrain. This state bifurcates in the positive u direction 
wherever ^'^ < 54, or equivalently 0q < 3g''/2. In the following we are interested in the modu- 
lational instability of this state. We suppose that this instability takes place at u = Uc and write 
u = Uc + 6^i>, where 5 ^ e is a new small parameter and z> = 0{1). In addition, we write 
A = Ao{l + A), Ox = V. Since the imaginary part of A decays to zero we take A to be real and 
write 

i = Ml + 6^A2 + 6^A-i + ..., V = 6Vi + 6^V2 + 6^V3 + .... (25) 
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Finally, since the localized states created at z/c are stationary, we set the time derivatives to zero 
and integrate Eq. (20) once obtaining 

V + ^\A\^ + C{S) = 0, (26) 

where C = Co + 5Ci + S'^C2 + S^Cs + . . . and the Cj are determined by the requirement that the 
average of V over the domain vanishes. 

Substitution of the above expressions yields a sequence of ordinary differential equations which 
we solve subject to periodic boundary conditions. At 0{1) we obtain Cq = —^\Ac\'^, where 
denotes the value of A^dXy = v^. At 0(5) we obtain the linear problem 

- 2z/,ii + 4iixx - = 0, \/i + 2e|A,|2ii = 0, Ci=0, (27) 

and conclude that A^ = An cos iX, Vi = Vn cos £X, provided that = -2f[{l - /(I - 

19^^/54)]. Note that this quantity is positive when 54 > > 54/19 (equivalently 3g^/2 > (pl > 
3g^/38) and that it vanishes in the limit £ — )■ 0, i.e., for infinite modulation length scale. 
At 0{5'^) we obtain 

- 2u,A2 + AA^xx - &2 = ^AV, + 3u,Al V2 + ^\Af{2A2 + Aj) + C2 = 0. (28) 
Thus A2 = A20 + A22 cos 2iX, V2 = V22 cos 2£X, where 

M2 = \a\,, V22 = -i\A,?Al,, (29) 
together with 

C2 = -\i\A,?Al,-2i\A,\^A2o. (30) 

Finally, at 0{5'^) we obtain 

- 2z/,i3 + ^Mxx - iV-i = 2vAi + i{AiV2 + isV^i) + //^(eiiia + A\), (31) 

and 

V, + i\A,\\2A, + 2A,A2) + C3 = - 3(1 _ (^2/54)) ^-^1- (^2) 

The direction of branching of solutions with A ^ follows from the solvability condition for ^3, 
i.e., the requirement that the inhomogeneous terms contain no terms proportional to cosiX. We 
obtain 

i> + 3fAl 
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It follows that the localized states bifurcate in the positive v direction (lower temperature) when 
54/19 < ^ < 54/13 and in the negative u direction when 54/13 < ^'^ < 54. The former 
requirement is equivalent to 3g^/38 < 0q < 3g'^/26, the latter to 3g^/26 < 0g < 3q^/2. These 
results agree with the corresponding results for an equation similar to Eq. (24) obtained by Elmer 
[79] (see also [22, 37]). 

In the above calculation we have fixed the parameter ^ (equivalently 0o) and treated r (equiv- 
alently u) as the bifurcation parameter. However, we can equally well do the opposite. Since 
the calculation is only valid in the neighborhood of the primary bifurcation at r = — 30q, i.e., in 
the neighborhood of ,^ = — V— 12r/g^, the destabilizing perturbation analogous to i/ > is now 
replaced by ^ > (cf. Fig. 1), implying that the direction of branching of localized states changes 
from subcritical to supercritical as decreases through = 54/13 (equivalently 0q decreases 
through 0Q = 3g'^/26 or r increases through r = — 9g''/26 ~ — 0.35g''), in good agreement with 
the numerically obtained value r ^ —0.41 (for q = 1) that is indicated in Fig. 7 by the horizontal 
line separating regions (iii) and (iv). 

We next consider the corresponding results in the case when the chemical potential fi is fixed. 
In this case the constants Cj vanish, and the C(5^) problem therefore has solutions of the form 
A2 = A20 + A22 cos 2iX, V2 = V20 + 1^22 cos 2eX, where 

A2o = -IaI„ V20 = ^\A,\^Al„ A22 = \aI„ V22 = -^\A,\^Al,. (34) 

The nonzero V20 contributes additional terms to the solvability condition at 0{5'^). The result 
corresponding to (33) is 



i> + 3f 



1 - (eV54) 



0. (35) 



1 - (19^754) _ 

This implies that the secondary LS branches are subcritical whenever the periodic branch is su- 
percritical (^^ < 54) and the secondary bifurcation is present (^^ > 54/19). These results are 
consistent with Fig. 13 and moreover predict that the LS states are absent for < 54/19, i.e., 
for 00 > -v/3/38g2 or r > -9g^/38 (cf. Fig. 13). This point is, of course, the thermodynamic 
tricritical point discussed in Sec. II B. 
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